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Abstract. In this paper we document the results of the 
study which led to the ripple correction and absolute cal- 
ibration algorithms applied to the high resolution spec- 
tra processed with the NEW SIPS software for the Final 
Archive of the lUE Project. In this analysis, based on a 
very large number of spectra, we find that both K and 
the a parameters (not only the former as previously be- 
lieved) vary with order number. This fact, together with 
the finding that the central peaks of the blaze function 
vary also as a function of the THDA temperature (for 
the SWF camera) and of the date of observations (for 
the LWP and LWR cameras), makes the ripple correction 
algorithm more complex than previously considered but, 
at the same time, considerably more reliable. As for the 
high resolution absolute calibration, the method followed 
is similar to the one implemented in lUESIPS. The in- 
ternal accuracy of the high resolution calibration is about 
4%. We note that the ripple correction and absolute cali- 
bration algorithms here described apply also to lUE data 
processed and distributed with the INES system. 

Key vi^ords: Methods:Data Analysis; TechniquesTmage 
Processing; Ultraviolet: General 



1. Introduction 

In this paper we address two fundamental aspects which 
influence, to different extent, the photometric quality of 
lUE high resolution spectra, i.e. the correction for the 
echelle blaze function (ripple correction), and the absolute 
calibration. The results of this study, which are based on 
about one thousand spectra, mainly of /f/i? standard stars. 
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have been implemented in the lUE Newly Extracted Spec- 
tra system {INES), a data archive and distribution facility 
developed by the lUE ESA project (see Wamsteker et al. 
1999), as well as in the NEWSIES data reduction pack- 
age (Nichols and Linsky 1996, Garhart et al. 1997). Full 
details on the INES system are provided by Rodriguez- 
Pascual et al. (1999, Paper I) for low resolution data, and 
by Gonzalez-Riestra et al. (1999a, Paper III) for high res- 
olution data. 

2. The NEWSIPS Ripple Correction 

The correction of lUE high resolution spectra for the 
echelle grating efficiency has for long been considered a 
critical area in the context of lUE data reduction pro- 
cedures. The INES system includes an upgraded ripple 
correction algorithm which is presented and discussed in 
the following. 

2.1. The method 

Let us indicate with Rm(A) the grating efficiency (blaze 
function) of the high resolution spectrographs, for a given 
order m and wavelength A. As shown by Ahmad (1981) 
and Ake (1981), this function is adequately represented 
by: 



Rm(A) 
where 



TT m a 



'x/x^ 



(1-Ae(m)/A) 



(1) 



(2) 



being Ac(m) the wavelength of the central maximum of the 
blaze function for order m, and a a parameter which can 
be shown to be, in the first approximation, inversely pro- 
portional to the Half Width at Half Maximum (HWHM) 
of the order considered: 



HWHM = 1.395 Ac(m)/(7r m a) 



(3) 
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The quantity 
K = mAc(m) 



(4) 



is called the 'ripple constant', but, as discussed in the fol- 
lowing, it turns out to be a function of order number for 
the lUE setup. 

Let us indicate with f(A) the observed net spectrum 
normalized to the exposure time. The shape of the blaze 
function Rni(A) can be determined from the observa- 
tions, knowing the inverse sensitivity curve of the camera 
(S~^(A)C(A), see Section^, and the absolute flux of the 
target F(A): 



R„,(A)=f(A)C(A)S-i(A)/F(A) 



(5) 



This equation implies that the blaze parameters can- 
not be derived from f(A) directly, because it is affected 
by the distortions introduced by the multiplying fac- 
tor C(A)S~-'^(A)/F(A), which is not only wavelength- 
dependent, but differs also from target to target. Not hav- 
ing taken this fact into account is one of the major causes 
for the inaccuracy of the ripple correction in previous re- 
leases of the lUE data reduction packages and, in particu- 
lar, for the suspected dependence of the ripple parameters 
on the spectral type of the source considered. 

Our approach to derive the ripple parameters is equiv- 
alent to using Eq. y, but is of much easier implementa- 
tion. It consists in normalizing the observed spectrum to 
its continuum. The continuum was determined, for each 
spectral order, in correspondence to the expected peaks 
of the blaze function (Grady and Garhart 1989), where 
i?in(A) ~ 1, and then interpolated over wavelength. 

Let us indicate with Y(A) the net spectrum extracted 
from order m after normalization to the continuum, and 
with T(m,A) its analytical representation. We can write: 



T(m,A) = A„,R^(A) 



(6) 



where Am is an adjustable parameter representing the net 
peak intensity of the blaze function (at A = Ac). Because 
of normalization to the continuum. Am ~ 1. The unknown 
blaze parameters a, Ac(m) and Am were determined itera- 
tively, for each order, via a non-linear least-squares fitting 
algorithm, i.e. by minimizing X] (^i — Yi)^. All the data 
points were given the same weight except those affected 
by saturation, ITF extrapolation, reseau marks or parti- 
cle events, which were discarded. Spectral absorption fea- 
tures were also discarded through an automatic rejection 
procedure. An example of the fitting for the LWP05471 
spectrum of the standard star BD+28 4211 is shown in 
Fig. 0. 

We stress that the ripple parameters should be derived, 
as done here, in vacuum wavelengths and in the velocity 
scale of the lUE spacecraft. 
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Fig. 1. Example of fitting Eq. || to a portion of a LWP 
spectrum of BD+28 4211. 



2.2. The data base 

For this study we have selected 516 good quality high res- 
olution spectra of the lUE calibration standards obtained 
between 1978 and 1995, and namely: 182 spectra for the 
SWP camera, 296 for the LWP and 138 for the LWR. 
This represents the near totality of data of this kind avail- 
able. As input data we used the MHXI net spectra as 
extracted with NEWSIPS. Among the most frequently 
observed standards we quote BD-l-28 4211, BD-l-75 325, 
HD60753, HD3360, G191 B2B, NGC246, HD120315 and 
HD93521. In addition, to test specific aspects of the ripple 
correction algorithm, it was necessary to introduce about 
200 net spectra (also extracted with NEWSIPS) of targets 
having heterogeneous spectral types. The large majority 
of the spectra were obtained through the large entrance 
apertures of the lUE spectrographs. A few small aper- 
ture spectra were also considered, whenever available, to 
cross-check the applicability of the large aperture ripple 
correction algorithm to that case. For the long wavelength 
cameras, several overexposed spectra were also included to 
refine the ripple parameters at the shortest wavelengths, 
where the cameras sensitivity is low. 

2.2.1. SWP 

Dependence of K on order number: We find that the 
central positions of spectral orders Ac(m) are not sta- 
ble for SWP spectra, as commonly believed in the past 
but, for a given order, vary strongly as a function of the 
Camera Head Temperature (THDA). Given this, a con- 
venient way to proceed is to determine first Ac(m) from 
data obtained within a restricted range of THDA values. 
We used 52 spectra of the stars BD-|-28°4211 (13 spec- 
tra), BD-h75°325 (20 spectra), CD-38°10980 (5 spectra), 
HD60753 (12 spectra) and NGC246 (2 spectra) obtained 
the range 8.5 < THDA < 11.5°C. Being the mean value 
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of THDA for this sample, 10.04±0.70 °C, very close to the 
average operating conditions for the SWP camera during 
years 1978 to 1991, the corresponding central positions 
Ac(m) can be used to obtain a mean curve K=m Are/, 
where Xref = Ac(m,THDA=10.04) are the central wave- 
lengths to be taken as a reference for that temperature. 
A linear regression to the mean values for the four stars 
provides: 



K = mAic 



137508.316 + 2.44761341m 



(7) 



The associated uncertainty on the central wavelengths is 
0.24 A at 1400 A and 0.31 A at 1850 A. The ripple- 
corrected fluxes in the overlap region between adjacent 
orders are accurate to within 3% above 1700 A and 6% 
around 1200 A. However, the flux mismatch in the overlap 
region between adjacent orders caused by a systematic er- 
ror on the central wavelengths is twice the quoted errors, 
because if the flux at order m is overestimated, the flux 
in the overlap region of the adjacent order m-1 will be 
underestimated and viceversa. 

Dependence of Ac on THDA: To extend the validity of 
Eq. to any regime of THDA, we have used the whole 
set of input spectra and computed for each spectrum the 
wavelength difference between the observed central wave- 
lengths and those obtained from Eq. 0. A linear fit to 
the data over the range of orders 110-86, where the blaze 
function is narrower, and then better sampled all over its 
shape, provides: 



Ac ^ Arc 



-0.335775 -I- 0.0321729 THDA 



(8) 



with a standard deviation of 0.069 A. This equation 
quantifies the previously quoted dependence of the 
central positions on THDA. We have also looked for a 
possible dependence of the central wavelengths on the 
epoch of the observation, on the spectral type of the 
target and on the exposure level, but no correlation was 
found. 
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Fig. 2. Ripple parameters for the SWP camera. Form top 
to bottom: Variation of K with order number, dependence 
of AA on THDA, and dependence of a on order number 
and on focus step. 



Dependence of a on order number and focussing condi- 
tions: We find that the ripple parameter a varies linearly 
as a function of the order number. A linear fit to the data 
provides: 



a = 0.92628 - 7.890132 IQ-^m 



(9) 



with a standard deviation of 0.022. The corresponding un- 
certainty on fluxes in the regions midway between adjacent 
orders is 3.6%, irrespective of the order considered. Note 
that the errors in a do not produce any flux discontinuity 
in the overlap region between adjacent orders. 

No dependence of a on the date of observation or 
THDA was detected. On the contrary, there is a marginal 



indication of a slight decrease of a with increasing focus 
STEP parameter. This effect, if real, would imply that the 
width of the spectral orders becomes larger as STEP de- 
parts from the optimum focussing conditions. Since most 
lUE data are taken at optimum focussing conditions and 
the effect is any case marginal, it has not been imple- 
mented in the ripple correction algorithm. In Fig. we 
show the ripple constant K as a function of order number, 
the wavelength shift AA — Ac— Arcf as a function of THDA, 
and the a parameter as a function of order number and 
focus STEP. 

Given the strong dependence of the central wave- 
lengths of spectral orders on the THDA temperature, we 
have verified that the SWP ripple correction algorithm 
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Fig. 3. Energy distribution of BD+75 325 obtained from 
three SWP ripple corrected high resolution rebinned 
spectra obtained under very different THDA conditions. 
Fluxes are in units of erg cm"^ s~^ A~^. 



STEP were 10.18±1.81 °C and -2.84±0.85, respectively. 
Since the central positions were sufficiently stable in the 
quoted period of time, the values of K=mAc(m,T=1991.2) 
obtained from the individual spectra were averaged to- 
gether. A linear regression to these data provides: 



K = mArcf(TO) = 230868.17770 + 3.86260914m 



(10) 



where Are/=Ac(m,T=1991.2) are the reference central 
wavelengths for time T=1991.2. The standard deviation 
of the fit is 30.1, corresponding to an uncertainty on the 
central wavelengths of 0.30 A at 2300 A. The ripple cor- 
rected fluxes in a point midway in the overlap region be- 
tween adjacent orders are accurate to within 3.9% at 2200 
A, and 2.3% at 2800 A. 



Dependence of Ac on time: In this section we show that 
the central wavelengths Ac of a given order vary with time, 
being their deviations with respect to the mean values in 
1991.2 (see Eq. |l(]| ) a linear function of observing time. 
To define this dependence, we have selected 142 spectra 
obtained between 1981 and 1995, and derived, for each or- 
der, the wavelength difference with respect to the reference 
values for 1991.2. Since the shape of the K(m) function is 
linear at any date and with a very similar slope, the wave- 
length differences of individual spectra could be averaged 
together over the orders 110 to 80, providing a mean wave- 
length shift A A = Ac -Aicf , which is plotted as a function 
of time in Fig. 0. A linear fit to these data provides: 



Ac - Arcf = 52.570796 - 0.0263910 T 
with a standard deviation of 0.10 A. 



(11) 



works well even when the THDA values depart consid- 
erably from the average operating conditions. As an ex- 
ample, we show in Fig. S, three spectra of BD+75 325 
obtained with very different THDA values (6.5, 9.5 and 
14.2): no periodic flux fluctuations (typical of bad ripple 
correction) are seen in the corrected spectra. 

2.2.2. LWP 

Dependence of K on order number: Contrary to the case 
of the SWP camera, the central wavelengths of LWP spec- 
tral orders Ac(m) are not sensitive to the camera head 
temperature THDA, but depend strongly on the date of 
the observations. 

It is then convenient to study first the dependence of K 
on order number for spectra obtained in a sufficiently short 
time interval. We have selected the period from year 1988 
to 1992 because many (86) good quality spectra of the lUE 
calibration standards were obtained. The mean observing 
date of this restricted sample was T = 1991.2 (fractional 
year), and the mean and rms values of THDA and focus 



Dependence of a on order number and focussing condi- 
tions: We find that the ripple parameter a shows a strong 
dependence on the order number. A linear fit to the data 
provides: 



a = 0.406835 + 0.01077191m - 5.945406 10 



-5m2 



(12) 



and a standard deviation of 0.00867. The corresponding 
uncertainty on fluxes in the overlap region midway be- 
tween adjacent orders is 1.4%. 

We have also investigated to which extent the cen- 
tral wavelengths depend on the camera head temperature 
THDA. To this purpose we have selected 194 spectra of 
calibration standards and computed the wavelength differ- 
ence AA=Ac — Are/ averaged over the order range 110-80. 
The results indicate that there is no correlation between 
AA and THDA, unlike the case of the SWP camera, as 
shown in Fig. ^. 

Similarly, no correlation has been found between the 
a parameter and the focussing conditions, nor between 
the blaze parameters and the exposure time or the energy 
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Fig. 5. Flux distribution of HD 3360 obtained from three 
LWP ripple corrected and flux calibrated high resolution 
rcbinned spectra obtained at very different epochs. Fluxes 
60 are in the same units as in Fig. 0. 
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year) and THDA, and the a parameter as a function of 
order number and focus STEP. 

We find that, in spite of the strong dependence of the 
central positions of spectral orders on the time of observa- 
tions (see Eq. 0), the quality of the ripple correction re- 
mains good even for spectra obtained several years apart, 
as shown in the example of Fig. pi. 



Fig. 4. Ripple parameters for the LWP camera. Form top 
to bottom: the K parameter as a function of order num- 
ber; the wavelength shift AA = Ac — Arcf as a function of 
observing date and THDA; the a parameter as a function 
of order number and focus STEP. Note the lack of corre- 
lation between AA and THDA, and between a and focus 
STEP. 



distribution of the star. The main results for the LWP 
camera are provided in Fig. 0, which shows the K param- 
eter as a function of order number, the wavelength shift 
AA=Ac — Aref as a function of observing date (fractional 



2.2.3. LWR 

Dependence of K on order number: As for the LWP cam- 
era, the central wavelengths of LWR spectral orders Ac(m) 
depend strongly on the date of the observations, but are 
insensitive to THDA variations. 

Following the same procedure used for LWP spectra, 
we have in first place determined a mean curve K(m)= 
mAc(m.) for a fixed epoch, to be taken as a reference. A 
restricted sample of 58 spectra of lUE calibration stan- 
dards obtained between 1982 and 1986 was used for this 
purpose. The mean observing date was T=1982.6 and the 
mean values of THDA and focus STEP were 13.40±1.70 
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Fig. 6. The K and a parameters for the LWR camera as 
a function of order number. The least-square polynomial 
fit to the data as from Eq. 13 and Eq. |l^, respectively, are 
also shown. 

°C and -1.3±0.9, respectively. The central wavelengths of 
spectral orders were then averaged together to obtain a 
mean value of K=mAc(m,T=1982.6). A fourth order poly- 
nomial fit to the results provides: 

K = A + B m + C m^ + D m'VEm'' (13) 

where 

A=0.281749635 10^ 
B=-0. 223565585 10^ 
C=0.365319482 10^ 
D=-0. 262477775 
E=0.701464055 lO^^ 

The standard deviation of the fit (18.06) corresponds to 
an uncertainty on the central wavelengths of 0.18 A at 
order 100. 

We find that the central wavelengths of LWR spectral 
orders Ac(m) do not show any dependence on the camera 
head temperature THDA, as for the LWP camera. 

Dependence of Ac on time: The central wavelengths ob- 
tained through Eq. 13 are strictly valid only for the 
T=1982.6 reference date. To study the dependence of the 
central wavelengths on time we have used 86 LWR spec- 
tra of lUE calibration standards obtained between 1978 
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Fig. 7. Examples of how the wavelength shift AA=Ac — 
Are/ for the LWR camera depends on the time of obser- 
vations for four different ranges of orders. The coefficients 
of the least-square polynomial fit to the data are given in 
Table 0. 



and 1994, and namely: 17 spectra of HD3360, 12 spec- 
tra of HD34816 and 19 spectra of HD120315. The results 
indicate that, differently from the LWP camera, the de- 
pendence on observing time is generally non-linear, and it 
varies across the camera faceplate. We have computed the 
mean value AA(m)=Ac(m)-Arcf (m) in 10 overlapping win- 
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Table 1. Coefficients to determine the wavelength shifts 
for the LWR camera 



HD 3360 



Order 


a(mi) 


b(mO 


c(mj) 


115 


35736.699 


-35.970790 


0.0090515542 


111 


62433.199 


-62.859885 


0.0158223297 


107 


53287.133 


-53.638642 


0.0134980459 


103 


42742.709 


-43.014583 


0.0108219635 


99 


28040.843 


-28.206111 


0.0070930274 


95 


10463.439 


-10.501169 


0.0026347077 


91 


6223.1919 


-6.2320255 


0.0015601476 


87 


4478.2512 


-4.4662128 


0.0011134034 


83 


112.32014 


-0.00566704 





79 


156.16074 


-0.0787367 






dows, each containing seven consecutive orders centered 
at 



mi = 115-4(i-l) with i= 1,10 



(14) 



The averages were made from order mj-S to order mj-f-S. 
A quadratic fit to the data is appropriate in all cases ex- 
cept near the long wavelength end of the camera, where 
a linear behaviour is found. In summary, the wavelength 
shift AA = Xrmc — Arcf Can be represented as: 



AA(mi, T) = a(mi) + b(mi)T + c(mi)T2 



(15) 



where T is the date of observation (fractional year). The 
coefficients applicable to individual orders are obtained 
by linear interpolation of the data in Table |^. For orders 
greater than 115 and smaller than 79 the coefficients for 
m — 115 and m = 79 should be used, respectively. 

Dependence of a on order number and focussing condi- 
tions: As in the case of the LWP camera, we find that 
the ripple parameter a has a strong dependence on the 
order number. A polynomial fit to the data provides: 



1.360633 - 4.252626 lO^^m 
(m= 119- 101) 
3.757863 - 0.0640201 m 

(m ^ 100 - 74) 



(16) 



3.5664390 lO^'^m^ 



with rms errors of 0.02 and 0.01, respectively. 

No significant correlation was found between the a pa- 
rameter and the focussing conditions. 

In Fig. P we show K and a as a function of order num- 
ber for the mean date 1991.2, and in Fig. H the mean 
wavelength shift AA = Ac — Arcf as a function of observing 
date. 

It is interesting to note that, in spite of the strong 
time dependence of the central positions of spectral orders 
(see Eq. |l^), the LWR ripple correction algorithm provides 
good results even if applied to spectra obtained several 
years apart, as shown in Fig. |8[ 
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Fig. 8. Flux distribution of HD 3360 obtained from three 
LWR ripple corrected and flux calibrated high resolution 
rebinned spectra obtained several years apart. Fluxes are 
in the same units as in Fig. g. 

2.2.4. The ripple correction algorithm 

In conclusion, the blaze function for lUE high resolution 
spectra processed with NEWSIPS remains analytically 
defined by Eq. n^, where the a parameter is given in Eqs. 
|, |l| and ^ for SWP, LWP and LWR, respectively, and 
the central wavelengths are computed according to: 

SWP: 
Ac(to) = 137508.316/m + 0.0321729THDA -1-2.111841(17) 
which is the result of combining Eqs. and g. 

LWP: 

Ac(m) = 230868.1770/m - 0.0263910 T + 56.433405 (18) 
which is the result of combining Eqs. HO^ and O. 

LWR: 

Ac(to) = K(m)/m + AA(m, T) (19) 
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where K is given by Eq. |T3| and AA(m,T) is computed 
according to Eq. |l^. 

The ripple correction algorithms were derived using 
spectral orders 119 to 72 for the SWP camera, 118 to 74 
for the LWP camera and 118 to 76 for the LWR camera. 
Extrapolation of the algorithms to higher and lower orders 
provides in general satisfactory results. 

3. The absolute calibration of NEWSIPS high 
resolution spectra 

For the absolute flux calibration of NEWSIPS high reso- 
lution spectra we have followed the method described in 
Cassatella et al. (1994). Let us indicate with N(A) the 
ripple-corrected high resolution Flux Numbers, normal- 
ized to the exposure time. The corresponding absolute flux 
can be determined from 



F(A) = C(A) S-^X) N(A) erg/cm^ s A 



(20) 



where S~^(A) is the low resolution inverse sensitivity func- 
tion appropriate to the camera considered, and C(A) is the 
so-called high resolution calibration function defined as 



C(A) = n(A)/N(A) 



(21) 



being n(A) the net Flux Numbers (i.e., not absolutely cal- 
ibrated), normalized to the exposure time, derived from 
low resolution observations of the same target. Because 
of the time-dependent sensitivity degradation of the cam- 
eras, the pairs of low-high resolution spectra used to de- 
termine C(A) should be obtained close enough in time. 
Alternatively, both n(A) and N(A) should be previously 
corrected for sensitivity degradation. This latter approach, 
here followed, has the advantage of increasing considerably 
the number of usable spectra. The correction for sensitiv- 
ity degradation of high resolution spectra has been made 
with the same algorithms used for low resolution spec- 
tra, as described by Garhart (1992, 1993) and Garhart et 
al. (1997). This procedure is justified in Paper III, which 
shows that high resolution spectra obtained even several 
years apart, once corrected in this way, provide very nearly 
the same flux repeatability performance as spectra ob- 
tained close in time. 

The high resolution spectra are first corrected for the 
blaze function, resampled in 2 A bins, and normalized to 
the exposure time in seconds. The spectra are then cor- 
rected for sensitivity degradation, for the THDA induced 
sensitivity variations, and for the camera rise time follow- 
ing the same algorithms used in the NEWSIPS process- 
ing of low resolution spectra (see Gonzalez-Riestra et al. 
1999b, Paper IV). Finally, the spectra of each target are 
averaged together to obtain a mean spectrum N(A). To ob- 
tain C(A), the mean high resolution spectra of a given tar- 
get are then divided by the mean low resolution spectrum 
of the same target (which are also corrected for tempera- 
ture effects, camera rise-time and sensitivity degradation). 



3.1. SWP 

To determine the high resolution calibration function C(A) 
we have used 28 SWP high resolution spectra of BD-l-28 
4211, 38 of BD-h75 325, 27 of HD60753, 6 of G191 B2B and 
13 spectra of CD-38 10980. We find that the repeatability 
errors on N(A) (after all the above mentioned corrections 
are applied) are typically 3-5%. These small errors confirm 
the validity of applying the low dispersion degradation 
rates to high dispersion spectra. 

The low resolution net fluxes n(A) of the above targets 
were obtained by averaging many low resolution spectra 
obtained during the 1990-1991 re-calibration period. The 
curve C(A) for the SWP camera is shown in Fig. ||. 

A third order polynomial fit to the data provides: 



C{X)^A + BX + CX'^ + DX^ 



where 



(22) 



A=1349.8538 

B=-2.0078566 

C=l. 10252585 lO'^ 

D=-2.0939327 lO"'^ 

with an standard deviation of 6.3. The above equation has 

been derived from data in the wavelength range 1175 to 

1950 A. 

The repeatability error on C(A) is 4%, irrespective of 
wavelength, which we take as the internal error of the high 
resolution calibration function. 

3.2. LWP 

To determine C(A) for the LWP camera we have used 25 
high resolution spectra of BD-h28 4211, 37 of BD-t-75 325, 
27 of HD60753 and 4 spectra of CD-38 10980. The re- 
peatability errors on N(A) reach the 4% level at 2400 A, 
but do not exceed 2-3% around 2800 A. Similarly to the 
case of the SWP camera, these small errors confirm the 
applicability of the low dispersion sensitivity degradation 
algorithm to high resolution data. The low resolution net 
fiuxes n(A) of the above targets were obtained by averag- 
ing many low resolution spectra obtained during the 1990- 
1991 re-calibration period, extracted with NEWSIPS and 
corrected for time-dependent sensitivity degradation ac- 
cording to Garhart (1993). The curve C(A) for the LWP 
camera is shown in Fig. n^. A linear fit to the measure- 
ments provides: 



C(A) = 251.383956 - 0.053935103 A 



(23) 



with a standard deviation of 3.49. The repeatability of the 
C(A) function is about 4%. The above equation has been 
derived from data in the wavelength range 1975 to 3150 
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Fig. 9. The high resolution calibration function C(A) for 
SWP spectra. Error bars are indicated. The thick line rep- 
resents the polynomial representation in Eq. 23 

3.3. LWR 

To determine C(A) we have used a total of 17 high and 23 
low resolution spectra of the calibration standards BD-l-28 
4211, BD+75 325, HD60753 and HD93521. 

We find that the repeatability errors on the net fluxes 
N(A) after the sensitivity degradation correction are of the 
same order as for the LWP camera, confirming once more 
the applicability of the low dispersion sensitivity degrada- 
tion algorithm to high resolution spectra. 

The resulting determinations of C(A) are reported in 
the bottom panel of Fig. |l^. It is interesting to note that 
the data points are well fitted by the same analytical rep- 
resentation as for the LWP camera (Eq. ^ ). The resid- 
uals correspond to an rms error on C(A) of 5.3. A linear 
fit would provide about the same residuals (4.5), and the 
curve would only deviate from that of the LWP by 1.3%. 
The internal accuracy of the calibration function ranges 
from 5% below 2300 A to 3% at longer wavelengths. The 
wavelength range covered by the LWR high resolution ab- 
solute calibration is the same as for the LWP camera. 

Note that also Cassatella et al. (1994) found that the 
C(A) curve is the same for LWP and LWR data processed 
with lUESIPS. 



Fig. 10. The high resolution calibration function C(A) for 
LWP and LWR spectra. Error bars are indicated. The 
thick line represents the adopted linear representation of 
C(A) in Eq. ||, which is the same for the LWP and LWR 
cameras. 
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Fig. 11. The flux distribution of the lUE standard 
HD60753 obtained from the high resolution spectrum 
SWP 9093 (full line) is compared with the absolute fluxes 
of the same star (open boxes) . The high resolution spec- 
trum has been rebinned to the low resolution wavelength 
step. Fluxes are in the same units as in Fig. g. 



3.4. Examples of application t. 1 c i, 

nux scale (Paper IV j. Examples ot such a comparison are 

It is important to compare the fluxes obtained from high given in Figs. iri|, O and Ol referring to the stars IID60753, 

resolution spectra using the present method with the BD-l-28 4211 and BD-t-75 325 observed with the SWP, 

fluxes of the calibration standards which define the lUE LWP and LWR cameras, respectively. 
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Fig. 12. The flux distribution of the lUE standard 
BD+75 325 obtained from the high resolution spectrmir 
LWP22343 (full line) is compared with the absolute fluxes 
of the same star (open boxes). The high resolution spec- 
trum as been rebinned to the low resolution wavelength 
step. Fluxes are in the same units as in Fig. 0. 
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Fig. 13. The flux distribution of the lUE standard 
BD+28 4211 obtained from the high resolution spectrum 
LWR12975 (full line) is compared with the absolute fluxes 
of the same star (open boxes). The high resolution spec- 
trum has been rebinned to the low resolution wavelength 
step. Fluxes are in the same units as in Fig. 0. 

In the above examples, the agreement between high 
and low resolution fluxes is within the 4% repeatability 
errors quoted above. 

It should be stressed that the present calibration is 
applicable to both continuum and emission line sources. 
This is confirmed by line emission measurements in sev- 
eral pairs of low-high resolution spectra of emission line 



sources. As an example, we show in Fig. 14 a low and a 



high resolution spectrum of the recurrent nova RS Oph 
taken very close in time. 
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1600 
Wavelength 



Fig. 14. Comparison between the flux calibrated high and 
low resolution spectra of the recurrent nova RS Oph. Dia- 
monds mark flagged pixels (saturated or affected by reseau 
marks) in the low resolution spectrum. Fluxes are in units 



of 10 ^^ erg cm ^ s 



Another test made was to verify the accuracy of the 
absolute calibration in the overlap region around 1950 A 
between the SWP and the LWP and LWR cameras. We 
find that, in this region, the short and long wavelength 
cameras agree to within 2 to 6% on average. As an exam- 
ple, we show in Fig. Ua the overlap region for two pairs of 
spectra of HD60753. 

The good match between short and long wavelength 
high resolution spectra can also be deduced from Fig. y^ 
which shows a combined SWP-LWP spectrum of the Wolf 
Rayet star HD152270. 

We note that residual non-linearity effects in the In- 
tensity Transfer Function, especially in the case of under- 
exposed spectra, or spectra near the saturation limit, can 
occasionally cause a fiux mismatch between the short and 
long wavelength cameras. 

4. Conclusions 

We have shown that a reliable correction for the echelle 
blaze function can be obtained for lUE high resolution 
spectra processed with NEWSIPS. We find that not only 
the K parameter, but also the a parameter varies with or- 
der number, and that the central peaks of the blaze func- 
tion deviate from the reference positions in a way which 
depends on the value of the THDA temperature for the 
SWP camera, and on the observation date for the LWP 
and LWR cameras. Having taken into account all these 
dependences (and not only the dependence of K on order 
number as in previous investigations), allows the present 
algorithm to successfully correct most spectra. According 
to our experience, the algorithm might fail in a few cases 
(probably less then 5-10%). The identified cases of failure 
correspond to images in which the interorder background 
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Fig. 15. The combined SWP-LWP and SWP-LWR spec- 
tra of HD60753 are shown in the overlap region around 
1950 A. The original data have been rcbinned in 0.5 A 
steps. Fluxes are in the same units as in Fig. |^. 




level is not properly evaluated with the NEWSIPS auto- 
matic extraction procedure or to cases in which the target 
is not acquired (purposely or not) in the center of the large 
entrance apertures of the lUE spectrographs. In this re- 
spect, we note that a displacement of 1 arcsec along the 
dispersion line corresponds to a velocity shift of the wave- 
length scale of about 5 km/s, irrespective of the camera 
used. For example, a displacement of 5 arcsec, or 0.125 A 
at 1500 A would produce a flux mismatch in the overlap 
region midway between adjacent orders of about 4%. 

As for the absolute calibration of lUE high resolution 
spectra, we have shown that its internal accuracy is about 
4% and that it applies both to continuum and emission 
line sources. The flatter shape of the C(A) function found 
here is due to the improved background extraction proce- 
dure in NEWSIPS (Smith 1999). 
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Fig. 16. Combined SWP-LWP spectrum of the Wolf 
Rayet star HD152270 after application of the ripple correc- 
tion and absolute calibration algorithms here described. 
The high resolution spectra have been resampled to 
the low resolution wavelength step for display purposes. 
Fluxes are in units of 10^^^ erg cm^^ s^^ A^^. 



